Field Dependence of the Ground State in the Exotic Superconductor CeCoIn 5 : 
a Nuclear Magnetic Resonance Investigation 
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We report 115 In nuclear magnetic resonance (NMR) measurements in CeCoIns at low temperature 
(T « 70 mK) as a function of magnetic field (-Ho) from 2 T to 13.5 T applied perpendicular to the 
c-axis. NMR line shift reveals that below 10 T the spin susceptibility increases as \fH~o. We associate 
this with an increase of the density of states due to the Zeeman and Doppler-shifted quasiparticles 
extended outside the vortex cores in a d-wave superconductor. Above 10 T a new superconducting 
state is stabilized, possibly the modulated phase predicted by Fulde, Ferrell, Larkin and Ovchinnikov 
(FFLO). This phase is clearly identified by a strong and linear increase of the NMR shift with the 
field, before a jump at the first order transition to the normal state. 

PACS numbers: 74.70.Tx, 76.60.Cq, 74.25. Dw, 71.27. +a 
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Applying a magnetic field to a superconductor is a 
powerful way of revealing the complexity of this macro- 
scopic quantum state of matter. Even though the ef- 
fects of the field in conventional type-II superconductors 
are well established, in a wide range of systems (such as 
the high-T c , heavy fermions, and organic superconduc- 
tors) the consequences are far from being understood. 
One possibility is that the magnetic field may affect 
the competition between superconductivity and antifcr- 
romagnetism (or other competing orders), as suggested 
by several experiments and theories in high-T c super- 
conductors 0,0. Particularly in the case of a d-wave 
superconductor, study of the low energy excitations in 
the vortex state can provide useful information about 
the nature of the underlying ground states. The local 
density of states (DOS) of these excitations can be ef- 
fectively probed by measuring the magnetic field depen- 
dence of the NMR shift as presented here. Single crystals 
of heavy-fermion superconductor CeCoIn 5 are an ideal 
material for such a study owing to their high purity, most 
likely d-wave gap symmetry 3, 4] , and entirely accessible 
magnetic field - temperature (T) phase diagram. 

Furthermore, an applied magnetic field (Hq) may in- 
duce a novel superconducting (SC) state in which the 
momentum of the Cooper pairs is not equal to zero, but 
becomes finite and proportional to Hq, as predicted by 
Fulde, Ferrell, Larkin, and Ovchinnikov (FFLO) In 
this state the SC order parameter oscillates in real space. 
The FFLO phase is expected to occur in the vicinity 
of H C 2 when Pauli pair breaking dominates over orbital 
effects. The search for this exotic SC phase has at- 
tracted much of current interest @, 0, H, 0, and there 
is experimental evidence that it is realized in CeCoIns 
S Ell EH- Moreover, CeCoIn 5 is close to an antifer- 
romagnetic (AF) instability, offering thus a unique pos- 
sibility to study the competition between FFLO and AF 
instabilities. It was recently shown that a substitution of 
10 % of Cd for In could generate AF droplets 12|. Even 



in pure CeCoIns, NMR has revealed the existence of some 
form of magnetism within the possible FFLO phase [l3| . 

Here we report detailed 115 In NMR measurements as 
a function of H from 2 T to 13.5 T at T = 73 mK. We 
establish that the shape of the NMR spectra is deter- 
mined by the inhomogeneous field distribution of a vor- 
tex lattice (VL) for Hq < 10 T. For fields above 10 T, 
where a possible FFLO state (henceforth referred to as 
"modulated" or mSC state) is stabilized, the NMR spec- 
tra broaden well beyond what is expected on the basis 
of the VL distribution. Magnetic field dependence of 
the shift reveals explicit effects of the Zeeman energy on 
quasiparticles in this Pauli-limited superconductor. Be- 
low 10 T, the spin susceptibility increases as y/Ho. This 
increase is related to the excess DOS due to the Zeeman 
and Doppler-shifted nodal quasiparticles, confirming the 
presence of gap nodes in the SC state of CeCoIn 5 . Above 
10 T, the spin susceptibility is strongly enhanced and in- 
creases linearly with Ho, as expected in the FFLO state 
[3| . If instead this enhancement is ascribed to field in- 
duced magnetism, our data place a lower-bound on the 
spatial extent of the magnetic regions. 

We have used high quality single crystals of CeCoIns 
grown by a flux method [15| . Here, the spectra of the ax- 
ially symmetric In(l) site in H || [100] (aligned to better 
than 2°) are reported. They were recorded using a cus- 
tom built NMR spectrometer and obtained, at each given 
value of Hq, from the sum of spin-echo Fourier transforms 
recorded at constant frequency intervals. The magnetic 
shift was determined by the diagonalization of the full 
nuclear spin Hamiltonian and by subtracting the orbital 
contribution of ~ 0.13 %. The RF coil was mounted into 
the mixing chamber of a 3 He/ 4 He dilution refrigerator 
while a variable tuning delay line and matching capaci- 
tor were mounted outside the NMR probe to allow for a 
wide frequency/field coverage. The 63 Cu NMR on cop- 
per nuclei from the RF coil was used to determine the 
exact value of Hq. The sample was both zero-field and 
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field-cooled and no visible influence of the sample's cool- 
ing history on the NMR spectra was detected. In order 
to avoid heating of the sample by RF pulses we used very 
weak RF excitation power [161 ] and repetition times of the 
order of ten seconds. 

In the Abrikosov SC state (henceforth referred to as 
the "uniform" uSC state), at low T, vortices tend to 
form a solid lattice [l7| resulting in a spatial distribu- 
tion of magnetic fields. This distribution is reflected in 
the NMR lineshape, as shown in Fig. [lji), and can be cal- 
culated by solving the Ginzburg-Landau (GL) equations 
[13, EH. We calculated the spectra employing Brandt's 
method [HI, valid at any value of Hq between the lower 
and upper critical fields. The input parameters for the 
calculation are the Hq, coherence length (£), penetration 
depth (A), and VL geometry (the angle a). For a given 
value of Hq, a \ 2 minimization of the difference between 
the calculated and measured spectra, normalized to their 
corresponding areas, was performed with £, A, and a as 
variational parameters. As the result is not very sensi- 
tive to the variation of a parameter, its value was set 
to a = 90°, i.e. a square vortex lattice. At Hq = 4T 
the minimum \ 2 was achieved for £ = 34 ± 10 A and 
A = 1580 ± 120 A, values consistent with earlier reports 
[l9l |. The spectra were then simulated for all field val- 
ues 2T < Hq < 11.6 T using this single set of param- 
eters. In order to make a quantitative comparison, the 
full width at half maximum (FWHM) and the square 
root of the second moment (Vo^) of the simulated and 
measured spectra are plotted in Fig. [Tb) as a function of 
Hq. The agreement between the simulation and the ex- 
periment in the uSC state for Hq < 10 T is excellent. To 
our knowledge this is the first time that the comparison is 
performed over such a wide field range. Within the sen- 
sitivity of our measurements, these results indicate that 
for H | I [100], ranging from 2 T to 10 T, the GL model 
provides an adequate description of the vortex state even 
for this Pauli-limitcd superconductor, contrary to the re- 
sults of neutron scattering study for Hq perpendicular 
to the planes [2(|. This discrepancy might stem, besides 
Hq orientation, from the fact that the probability distri- 
bution, i.e. NMR spectrum, is insensitive to the Pauli 
effects which are enhanced near the vortex cores (2lj . 

For fields above 10 T, the measured spectra broaden 
beyond what is expected from the calculated VL distri- 
bution in the GL model. We have found that it is impos- 
sible to account for the high field linewidth broadening 
regardless of the £ and A values used. An appreciable 
discrepancy begins at Hq = H* that corresponds to the 
transition field from the uSC to the mSC state @, fl6| . 
Therefore, we demonstrate that at low T the transition to 
the mSC state is identified by the onset of the additional 
asymmetric NMR linewidth broadening. This implies the 
presence of an additional modulation of the internal field 
in the mSC state. The lineshape contains a single peak, 
which is consistent with an incommensurate modulation 
of the internal field along two spatial directions 1231 . as 
expected for an FFLO state in a <i-wave SC [2^, |24| . 



a) Frequency [MHz] b) 




3.98 3.99 4.00 4.01 2 4 6 8 10 12 



Hint m H m 

FIG. 1: (Color online) a) In(l), (-1/2 <-> -3/2) transition, 
spectra (red) at T — 73 mK and ffo = 4T compared to cal- 
culated VL lineshape (blue), as described in the text, as a 
function of the internal magnetic field (Hi nt ). The arrows 
indicate Hint corresponding to the saddle point (H s di) and 
Ho fields, respectively, b) FWHM and \/a^ of the mea- 
sured In(l) spectra (filled symbols) as a function of Hq. 
Solid lines are the calculated values obtained as described 
in the text. The dashed line denotes the field H* . 

We now proceed to the discussion of the main result, 
the Hq dependence of the low T shift data illustrated 
in Fig. [2] Two phase transitions are clearly discerned. 
The first order phase transition from the normal to the 
mSC state at Hq = H C 2 is characterized by the sizable 
discontinuity of the shift. The continuous second order 
phase transition from the uSC to mSC state occurs at 
Hq = H* . Presented data provide the first clear NMR 
signature of this phase transition [l3|, [l6j]. In both SC 
states, the shift increases with increasing Hq. However, 
above H* the field dependence of the shift is significantly 
enhanced. In the following, we will first address variation 
of the shift in the uSC. 

In a SC with the c?-wave gap symmetry in the excita- 
tion spectrum, at low T the quasiparticle excitations are 
restricted to nodal regions [25j. These low energy quasi- 
particles can be probed with NMR shift measurements. 
Specifically, the shift is proportional to the DOS (N(Ej) 
averaged over energy E in a r ang e of the order of fcgT 
around the Fermi energy, Ef [26(. Thus, at sufficiently 
low T the shift is proportional to the DOS at the Fermi 
level, i.e K ~ (N(E)) = Np. This implies that in the 
SC state only the DOS in the regions around the nodes 
of the gap contributes to K. Near the nodes the DOS 
depends linearly on quasiparticle excitation energy, for 
energies less than half of the gap magnitude. Therefore, 
the low T shift is given by, 

K~{N F )<x{\e + E z +E D \), (1) 

where (...) denote the average over four nodes in the k- 
space, and e ps fc^T, Ez = ^ef/Ho, and Ep, = vp ■ p s 
are the thermal, Zeeman, and Doppler energy terms, re- 
spectively. The Doppler term originates from the shift 
of the excitation energies of the nodal quasiparticlcs, 
with Fermi velocity Vf , moving in the superflow field 
with momentum p s . Both Ez and Ep> depend on Hq. 
The Hq dependence of the local p s , and thus Ep, at any 



3 




0.8 r 

' i i i i i i 

2 4 6 8 10 12 14 

H [T] 

FIG. 2: (Color online) In(l) magnetic shift (open circles) 
as a function of Ho at T — 73 mK. The solid circles are the 
shift corrected for diamagnetic contribution, i.e. shift with 
respect to the average magnetic field in the VL unit cell. 
The dotted lines are guide to the eye. The darker (red) 
solid line is the calculated shift, as described in the text. 
The lighter solid lines are a fit of the data, for Ho < 10 T, 
to Eq. [2] The dashed line is a linear fit to the data for 
Ho > 10 T. Inset: Magnetic field dependence of the local 
In(l) shift Ki oc at the saddle point field. The solid line is 
calculated shift, K = 1.04%+77.1- {e+E z +E Dsdl ) VoeV' 1 . 

point in the real space of the VL, can be calculated ex- 
ploiting Brandt's algorithm fl8| . The same algorithm 
is then used to directly calculate the H dependence 
of the shift averaged over the real space of a VL unit 
cell from A = Kq + K e ({\e + E z + E D \)). The follow- 
ing input parameters were used: vpH = 49.1 meVA [27| . 
E z = 0.5fi B H , K Q = 0.97%, and A e = 85.24 %eV"\ 
determined so that A at H°Z, h , the orbital H c2 = 38.6 T, 
equals the normal state shift. Excellent agreement with 
the data for Hq < 10 T is obtained with this simple 
model, as illustrated in Fig. [21 

We point out that the shift has been determined from 
the frequency corresponding to the first moment, that is 
the average frequency, of the spectra and thus reflects the 
real space average of K over a VL unit cell. In Fig. [2] we 
have thus presented both the shift data with respect to 
the applied field H , as well as the same data calculated 
with respect to the field averaged over the VL unit cell 
as obtained from the Brandt's algorithm. In our calcula- 
tion of the average A in the uSC state, we neglected the 
Pauli paramagnetic effects and their variation across the 
VL. This is justified by the fact that these effects are in- 
significant outside and are enhanced only near the vortex 
cores. Thus, mainly due to the large volume contribution 
from the outside of the cores, the average DOS is not no- 
tably affected by the Pauli paramagnetism [2l[ . The field 
dependence of the data is fully explained by the increase 
in the dominant Zeeman and average Doppler energy of 
the quasiparticles. Note that this behavior is in sharp 
contrast with that in a conventional SC without nodes, 



where the Doppler term has a negligible effect. Thus, our 
results indicate that CeCoIns is a SC with nodes in the 
gap, which is most-likely of d-wave (d x 2_ y 2) symmetry. 

The observed field dependence is also consistent with 
Volovik's prediction [25j for a <i-wave SC that the average 
DOS cx yfHu/H c2 . As shown in Fig. H H < 10 T data 
is well fitted to 

K = K + 13 K n ^H /H°* b , (2) 

where is fixed to 38.6 T and K n = 1.77% is the 

normal state shift. The fit parameters are Kq — 0.90 ± 
0.01 % and (3 = 0.33 ± 0.02. In the limit H -> 0, K 
remains finite, at nearly half of the normal state shift. 
The finite shift, that is the field independent constant 
contribution Ao, may be attributed to the multi-band 
nature of SC in CeCoIns. The contribution Kq would 
then originate from normal quasiparticles in the small 
gap band [2&| . 

In order to verify the significance of the contribution 
of Ed to the quasiparticle energy, we have also extracted 
the local shift from the peak of the spectra which corre- 
sponds to the nuclei positioned at the saddle point of the 
field distribution, that is at the point in real space mid- 
way between two vortices. At this point where H = H sc u, 
the local p s , and thus Ed, increases very slightly with in- 
creasing Ho (for Hq <§; H° 2 b ) since it is not influenced by 
geometry effects such as the change in the vortex num- 
ber with varying Hq. Thus, the local shift at the sadle 
point should exhibit weaker field dependence than the 
shift averaged over the entire VL. This is indeed the case 
as shown in the inset to Fig. [21 where the local shift 
at H s di is displayed. This shift corresponds to that of 
the peak of the spectrum calculated with respect to the 
H s di, obtained from Brandt's algorithm. Here the ob- 
served field dependence arises solely from the increase of 
the quasiparticles Zccman energy. 

We next consider the field dependence of the shift in 
10T < H o < H r2 . As apparent in Fig. [21 in this regime 
K increases linearly with Ho- Contrary to Q, no evi- 
dence of discontinuous jump in A, indicating the transi- 
tions between different Landau levels within the FFLO 
state, is observed. Our data can be well described by 

A = A h + h K n (H /H c2 ) , (3) 

where H c2 is fixed to 11.7 T, A„ = 1.77% is the normal 
state shift, h = 0.77 ± 0.02, and A# = 0.03 ± 0.01. The 
latter parameter shows that as Hq — + 0, A extrapolates 
to zero. The rate of increase of the shift (AA/AiJo) is 
five times higher than that in the uSC for Hq < H* . Af- 
ter careful consideration of all the possibilities regarding 
the values of Aq and the Doppler term, we conclude that 
this large increase can be ascribed to an enhancement of 
Ez and its dominance over Ed in the entire VL unit cell. 
For several reasons the fulfillment of Ez > Ed condition 
is very likely in the FFLO state. First, the internal field, 
due to paramagnetic effects, can be large enough in the 
FFLO state to exceed the Ed term everywhere in the VL 
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unit cell. Second, it is possible that Ed is suppressed by 
either the opening of a subdominant gap, or a change in 
the gap structure so that vp and ps become nearly or- 
thogonal. The gap opening scenario is very unlikely since 
K increases in the FFLO state contrary to expectations 
in the presence of a gap [13, 0|. On the other hand, 
changes in the gap structure are expected in the FFLO 
state. The observed AKf AHq could reflect the fact that 
the DOS is no longer simply proportional to E. That 
is, the DOS dispersion relation changes and/or an extra 
structure, such as bound states, is induced by Hq. This 
is consistent with the FFLO scenario [23[ but further cal- 
culations are required for a quantitative comparison. 

Nonetheless, a magnetic origin of the mSC phase or 
coexistence of some mag netic order with the FFLO state 
cannot be excluded [111 . 1 1 31 ] . Given the a priori antago- 
nistic relationship between SC and magnetism, it is likely 
that magnetism appears in the spatial regions where SC 
is suppressed as is the case in the vortex cores [2, IT^ . 
However, it was shown that the existence of the local 
moment magnetism requires pair coherence 11]. The 
absence of magnetism in the normal state could be at- 
tributed to the Kondo screening, acting on a short length 
scale, of Ce local moments. In the SC state the Kondo 
screening becomes ineffective, since quasiparticles have 
very small momentum inhibiting short range screening, 
giving rise to magnetism. Besides, a magnetic phase 
could be stabilized only in high fields when there is a 
sufficient overlap between vortex cores, so that the cor- 
relations between magnetic regions can be established. 
In field of 10 T, the distance between neighboring vor- 
tices is ps 140 A. Thus, the long range magnetic order 
could be established in fields above 10 T if the radius 
of the magnetic cores is larger than 2£ ps 70 A. In this 
case, the field enhancement of the shift would be a con- 



sequence of the canting of local Ce moments. Varying 
H from 10.2 T to 11.7 T increases the shift from 1.22 % 
to 1.38 %. With a hyperfine coupling between In(l) nu- 
cleus and its 4 Ce neighboring atoms of A ps 1.2 T//x# 
(29| . the increase of the average local magnetic moment 
is Afi e ff = A(KH )/4:A ps 0.0075 Ms/Ce. This corre- 
sponds to 7.5 % of a typical normal state local Ce mo- 
ment (~ 0.1 hb), implying that for H* < Hq < H C 2, the 
increase of Hq induces a canting of ~ 7.5 % of the normal 
state local Ce moment. 

In conclusion, our essentially zero-T limit data provide 
the clearest NMR evidence of two phase transitions in 
the vicinity of H C 2 thus far. In the uSC phase, the NMR 
spectra can be nicely fitted to the calculated magnetic 
field distribution arising from the vortex supercurrents 
with SC coherence length (£ ss 34 A) and penetration 
depth (A pa 1580 A) as fitting parameters. To explain 
the observed increase of the spin susceptibility with Ho, 
we conclude that the dominant low energy excitations are 
the Zeeman and Doppler-shifted quasiparticles extending 
outside the vortex cores, implying that CeCoIn 5 is of d- 
wave gap symmetry. For the high-field (above 10.2 T) 
low-T phase we find that it cannot correspond to a simple 
VL rearrangement. It is consistent with an FFLO state 
with 2D incommensurate modulation of the quasiparticle 
density in which the spin susceptibility strongly increases 
as Hq- If the magnetism appears there in the spatially 
localized fashion, we show that magnetic cores extend on 
a length scale larger than 2£. 
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